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Diffusion thermopower of (Ga,Mn)As/GaAs tunnel junctions 

Ts. Naydenova, P. Diirrenfeld, K. Tavakoli, N. Pegard, L. Ebel, K. Pappert, K. Brunner, C. Gould, L.W. Molenkamp 

PhysikaUsches InsUtut (EP3), Umversitat Wiirzburg, 
Am Hubland, D-97074 Wiirzburg, Germany 
(Dated: October 31, 2011) 

We report the observation of tunneling anisotropic magnetothermopower, a voltage response to 
a temperature difference across an interface between a normal and a magnetic semiconductor. The 
resulting voltage is related to the energy derivative of the density of states in the magnetic material, 
and thus has a strongly anisotropic response to the direction of magnetization in the material. The 
effect will have relevance to the operation of semiconductor spintronic devices, and may indeed 
already play a role in correctly interpreting the details of some earlier spin injection studies. 
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The most pressing issues of modern information tech- 
nologies are power consumption and heat dissipation. As 
such, a proper understanding of the relationship between 
electrical and thermal effects in devices is essential to 
their design and operation. When these devices are mag- 
netic in nature, the study of this relationship is the topic 
of the field of spin caloritronics 1], which was strongly 
inspired by the initial reports of the spin Seebeck effect in 
metals Q and subsequently in magnetic semiconductors 
i. 

As a simple example of the importance of thermal ef- 
fects on device operation, consider the problem of the 
injection of spin polarized current into a nonmagnetic 
semiconductor. Spin injection is often measured using a 
four contact nonlocal technique 0, Hj, where current is 
passed between two adjacent contacts, and the chemical 
potential resulting from spin accumulation is detected at 
two further contacts which are situated next to, but not 
between the current contacts. Ideally, the signal from 
such measurements should produce a voltage at the de- 
tection leads, which is equal but opposite in sign depend- 
ing on the relative orientations of the contacts. In prac- 
tice this is almost never the case, as the measurements 
typically include a constant voltage offset, which it has 
become customary to neglect as an unimportant back- 
ground contribution 

This background voltage is a thermoelectric effect Q , 
which in most cases is isotropic to magnetic field and 
thus independent of the magnetization direction of the 
contacts. The subtraction of such a constant background 
is thus inconsequential to the experimental conclusions. 

As we show in this Letter however, in materials where 
strong spin-orbit coupling links the magnetic proper- 
ties to the density of states 01, a spin caloritronic phe- 
nomenon, which we dub tunneling anisotropic magne- 
totermopower (TAMT), appears. This not only may have 
interesting implications for device applications, but can 
also lead to misinterpretations if one fails to take it into 
consideration. Indeed, we show that thermal voltages 



detected in such materials, e.g. the ferromagnetic semi- 
conductor (Ga,Mn}As, closely resemble typical nonlocal 
spin- valve signals 8], even when no spin injector is used. 

The configuration we study looks at the thermopower 
across a tunnel junction between ferromagnetic and non- 
magnetic regions [9] . Specificly, our experiment examines 
diffusion thermopower through a low-temperature grown 
i-GaAs tunnel barrier between a magnetic (Ga,Mn)As 
layer at 4.2 K and a highly Si doped GaAs channel 
in which the electronic system is heated by an ac cur- 
rent. Previous reports of thermopower in (Ga,Mn)As 
[lo| focused on the much larger phonon-drag contribu- 
tion which dominates thermopower in the bulk of the 
material. By designing our structure with a p-n inter- 
face and applying the heat directly to the electron sys- 
tem [ll^lSj, phonon-drag terms do not contribute to the 
measured thermovoltages ■ We can thus extract 

the effects of magnetic anisotropics on the diffusion ther- 
mopower term. This carrier-related thermopower probes 
the energy derivative of the (Ga,Mn)As density of states 
(DOS) and thus shows a magnetic anisotropy similar 
to (Ga,Mn)As tunneling anisotropic magnetoresistance 
TAMR This is clearly distinct from the twofold 

anisotropic ma gne toresistance AMR-like symmetry ob- 
served in Ref. |1(| for the phonon-drag dominated See- 
beck coefficient. 

Our structure has a single magnetic (Ga,Mn)As con- 
tact (20 nm thick), which is separated from a highly 
doped Si:GaAs mesa (60 nm) by a 1 nm low-temperature 
GaAs tunnel barrier. The semiconductor layers are 
grown on a semi-insulating GaAs substrate and buffer 
by means of low-temperature molecular beam epitaxy. 
Subsequently, and without breaking the vacuum, the 
sample is transferred to a UHV electron-gun evapora- 
tion chamber where Ti/Au contacts are deposited on 
the (Ga,Mn)As layer. The sample is patterned into the 
shape depicted in Fig. [1] using a multistep optical lithog- 
raphy and etching process. Contacts to the GaAs channel 
are established by Ti/Au evaporation and lift-off. The 
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FIG. 1: Schematic of our TAMT device showing the geometry 
of all contacts as well as the crystal orientation. Heating cur- 
rent is passed through the Si-GaAs channel, and the thermo- 
voltage is detected between contacts 2 and 5. The remaining 
contacts are used to inject the heating current and monitor 
the weak localization. 



50 X 50 fim^ pillar is connected to a bonding pad through 
a gold air bridge fabricated using electron-beam lithog- 
raphy and a multiple acceleration voltage process [l3l • 

Transport measurements are performed at 4.2 K in a 
variable temperature magnetocryostat fitted with a vec- 
tor field magnet that allows for the application of a mag- 
netic field of up to 300 mT in any direction. As expected 
from our targeted electron heating method IJ, |l5| , the 
(Ga,Mn)As contact remains at 4.2 K throughout the ex- 
periment, as confirmed by verifying that the anisotropy 
fingerprints of the (Ga,Mn)As remains unchanged when 
applying a heating current, and comparing that to the 
known temperature dependence of such fingerprints [3| ■ 

The electron system in the GaAs channel is heated out 



of equilibrium with respect to the GaAs lattice |l4 Il5| 
by a 13 Hz ac current applied across the channel, with 
one of its contacts defining the system ground. The re- 
sulting temperature difference AT across the tunnel bar- 
rier generates a thermal voltage Vth [3| which can be 
detected between contacts 2 and 5, as shown in Fig. [1] 
The amplitude of this thermal voltage depends on AT, 
or equivalently on the electron temperature in the GaAs 
channel. Since we heat at sufficiently low frequencies that 
the electron temperature remains in quasistatic equilib- 
rium and follows the applied power, the thermal voltage 
will be related to the square of the applied voltage, and 
can be detected at the double frequency 2/ using stan- 
dard lock-in amplifier techniques. 

To establish a relationship between the temperature of 
the carrier system in the GaAs channel and the heating 
current, we record the magnetoresistance of the GaAs 
channel in response to an out-of-plane magnetic field for 
various applied voltages. The width of the weak local- 
ization feature (20| as a function of heating current is 
then compared to the width of the same feature taken 
at varying cryostat temperature using low current. Do- 
ing so, we find that a heating areal current density of 5.5 



A/m through the heating channel corresponds to an elec- 
tron temperature of ~ 22 K, to an estimated accuracy of 
about ±1.5 K. 




FIG. 2: Amplitude of the thermopower signal as a function 
of the temperature in the heating channel. A magnetic field 
of 300 mT is applied along the [100] crystal direction. 

Fig. [5] shows the measured thermopower C = 
for different thermal gradients AT recorded at a field 
of 300 mT along the [100] magnetic easy axis of the 
(Ga,Mn)As layer. Similar values are obtained for other 
magnetic field directions and amplitudes. The measured 
voltage is clearly temperature dependent and has a mag- 
nitude in the /iV range, comparable to the background 
voltages typically observed in nonlocal measurements. 

The intercept with the y-axis S (4.2 K) = hmAT->o C — 
^ yields a Seebeck coefficient of 0.5 /^VK~^. This num- 
ber represents the diffusion part of the Seebeck coefficient 
at 4.2 K and is 3 orders of magnitude smaller than the 
phonon-drag dominated Seebeck coefficient measured in 
(Ga,Mn)As by Pu et al. 

This value is consistent with an estimate of the dif- 
fusion Seebeck coefficient at 4.2 K from a simple model 
of detailed balance across a p-n interface under the as- 
sumption that the carrier distribution on each side of the 
interface is represented by a Fermi function at a differ- 
ent temperature 2l|. We further assume that the thin 



tunnel barrier at the interface in our device can be ap- 
proximated to be large compared the transport energy 
scales and that the tunnel probability is energy indepen- 
dent. We thus model it by a delta function. 

This leads to an expression for total current through 
the interface given by: 



Jtotal 



A 



Ev,(Ga,Mn)As —eVth 



F[E)dE 



(1) 



Ec.Ga 



where 



FiE) ^ DoaAs-ME) ■ D^Ga,Mn)AsiE " cVth) X 

ifcaAs-.S'iiE) - f(Ga,Mn)As{E - eVth)) ■ (2) 
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Here, DcaAs-.SiiE) and D(^Ga,Mn)As{E - eVth) are the 
GaAs and (Ga,Mn)As density of states, and fcaAs-.SiiE) 
and f(Ga,Mn)As{E — eVth) are the Fermi functions in the 
two regions. ^ is a constant which includes geometry 
terms as well as the tunneling probability, and which is 
of no importance to the analysis. Since the total current 
flowing must be zero at equilibrium, given the density 
of states of the two materials the above equation can be 
solved numerically for Vth- 

For an order of magnitude estimate we assume a simple 
parabolic band model. The carrier density in the GaAs 
layer is n = 4 • 10^*cm~^ (from Hall measurements) and 
we estimate a value oip = 1 • 10^"'^ cm"'^ for the (Ga,Mn)As 
(typical for layers with similar growth conditions). The 
effective mass in GaAs is known and for the (Ga,Mn)As, 
an average hole mass of O.Smg is used. This calculation 
leads to a value of 0.3 fiVK~^ for the diffusion Seebeck 
coefficient, consistent with the experimentally obtained 
value of 0.5 i^VR-^ 

This model also predicts the linear dependence of ( 
on temperature as expected in general for the diffusion 
Seebeck effect [i^j and is in qualitative agreement with 
Fig. [21 While a quantitative model would require a de- 
tailed treatment of the top four valence band levels, in- 
cluding their anisotropic effective mass, the simplified 
degenerate semiconductor model implying a Fermi level 
deep in the (Ga,Mn)As valence band is sufficient to cap- 
ture the essence of the results. 

While this linear dependence is expected from the dif- 
fusive contribution to thermopower, phonon-drag contri- 
butions are typically highly nonlinear and usually even 
nonmonotonic [l^. Our observation of linear behavior 
is therefore further evidence that we are measuring the 
intrinsic carrier thermopower. 




FIG. 3: Thermovoltage in a magnetic field of 300 mT as a 
function of field angle. 0° is along [010] crystal direction. 

Since the (Ga,Mn)As DOS is the only magnetic field 
dependent term in Eq. [51 thermopower measurements re- 
flect its response to a magnetic field. More specifically. 



the product of D(^Ga,Mn)As{E) with a Fermi function 
leads to a response which is proportional to the energy 
derivative of the DOS of (Ga,Mn)As at the Fermi level. 
This implies that the observed signal should refiect the 
properties of the magnetic anisotropy of the (Ga,Mn)As 
DOS. 

Figure [3l shows the measured thermovoltage in our 
sample versus the direction of a 300 mT magnetic field ro- 
tated in the sample plane and thus continuously perpen- 
dicular to the temperature gradient. The magnetic field 
angle is given with respect to the [010] crystal direction 
as denoted in Fig. [1] We observe four equivalent minima 
close to the [100] and [010] crystal axes and two sets of 
local maxima: a larger one along [110] and a smaller one 
along [110]. This symmetry is strongly reminiscent of 
(Ga,Mn)As TAMR measurements which map the sym- 
metry of the (Ga,Mn)As DOS [16|] and is consistent with 
the thermovoltage mapping the energy derivative of the 
(Ga,Mn)As DOS as it depends on magnetization direc- 
tion. 

Because it maps the derivative of the (Ga,Mn)As DOS, 
the symmetry of the diffusion thermopower term in Fig. [31 
is very distinct from the AMR-like symmetry [2^ [23| ob- 
served in Ref. [l3| for the phonon-drag dominated See- 
beck coefficient. AMR has a cos^O dependence on the 
angle 9 between current and magnetization and thus a 
twofold symmetry with a minimum in resistance when 
the magnetization points along the current direction and 
a maximum for magnetization perpendicular to the cur- 
rent. In our sample, the current in the GaAs heating 
channel is along the 90° direction, while the tunnel inter- 
face between GaAs and (Ga,Mn)As is not in the current 
path. The location of minima and maxima in Fig. [31 is 
crystal direction dependent and has no relation to the 
current direction, as confirmed on an identical structure 
with the heating channel oriented along [100], and thus 
rotated 90 degrees compared to Fig. [H 

For a more detailed investigation of the magnetic 
anisotropy, we record the thermopower during high res- 
olution magnetic field scans from -300 mT to +300 mT 
along many in-plane crystal directions (Fig. H]) . The elec- 
tron temperature of the GaAs heating channel is kept 
constant at 12 K during the whole measurement, while 
the lattice temperature in the whole sample, as well as 
the carrier population in the (Ga,Mn)As layer remain at 
4.2 K. The recorded thermovoltage is a function of the 
magnetization direction in the (Ga,Mn)As layer and can 
again be explained by considering the (Ga,Mn)As mag- 
netic anisotropy. 

When the magnetic field is swept from -300 mT to 
300 mT, the magnetization rotates towards an easy axis 
of the material, here towards 0° and 90° as the field ap- 
proaches mT. It reverses at small positive fields in a 
mostly two-jump switching event. As the magnetization 
rotates, the thermovoltage changes smoothly. The largest 
magnetization rotation angle and thus largest change in 
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FIG. 4: Thermovoltage as a function of magnetic field swept 
along various in-plane directions. Inset: Enlargement of the 
low field part. The 190° curve is offset by 0.15 fj,V to remove 
a small thermal drift. 
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FIG. 5: Polar plot of the magnetic field positions of the 
switching events Ifd and IIc2 in the thermopower signal. 



thermovoltage is observed when the magnetic field is 
swept along the 45° hard axes of the (Ga,Mn)As ma- 
terial ([110]). When sweeping the field along or close 
to a magnetic easy axis ( 0° or 90°), the magnetization 
stays along this easy axis direction, thus the thermovolt- 
age is almost constant throughout the sweep apart from 
the one- or two-jump switching event. The figure thus 
shows two distinct components of TAMT: a volatile con- 
tribution caused by the field holding the magnetization 
away from its easy axes at high fields and a nonvolatile 
contribution associated with the difference in density of 
states when the magnetization switches between the two 
easy axis. 

It is worth noting that the switching features in the 
thermovoltage curves in Fig. 2] could easily be mistaken 
for the detection of a nonlocal spin signal. As an example, 



the inset displays the 190° curve (field scanned from left 
to right) together with the 10° curve (field scanned from 
right to left) in spin- valve fashion showing a similar spin- 
valve signal as reported in Ref.jl]. 

We extract the first and second switching fields of the 
magnetization reversal, Hd and Hc2 respectively, and 
plot them radially as a function of the field orientation 
(fi in Fig. m This plot reveals the (Ga,Mn)As anisotropy 
components, just as in previous crystalline AMR ^25] and 
TAMR 16| measurements. The global easy axes corre- 
spond to the diagonals of the rectangular pattern of the 
first switching fields Hd. As typical for (Ga,Mn)As [26j . 
two small uniaxial anisotropy components along the [110] 
(elongation of the center square) and along the [100] (step 
features in the central shape) crystal direction are visible 
in addition to the main biaxial anisotropy. 

The behavior of the measured voltage is thus fully de- 
scribed by the magnetization direction in the (Ga,Mn)As 
layer, together with the thermovoltage dependence on the 
magnetization direction given in Fig. [3l 

In summary, diffusion thermopower can be measured 
over a tunnel barrier between a magnetic (Ga,Mn)As 
layer at 4.2 K and a heated carrier system in a GaAs 
mesa. This diffusive component of the Seebeck coef- 
ficient is in accordance with theoretical estimates, and 
much smaller than the phonon-drag dominated Seebeck 
coefficient reported by Pu et al. 10]. Its magnetic field 
dependence is given by the derivative of the (Ga,Mn)As 
DOS and is also clearly different from the AMR-like sym- 
metry observed in Ref. ,10]. Moreover, this dependence 
of the thermopower on derivative of DOS could prove a 
useful tool in mapping out the details of the DOS, espe- 
cially in samples near the metal-insulator transition. 
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Note added - During the review of this Letter, three 
related papers (27l - [29j | appeared on spin thermopower 
across tunnel junction. These relate to metallic systems, 
and do not show anisotropic behavior. 
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